% g = f_g_new(bfu(k).data);
% J = f_dgdu_new(bfu(k).data);

g = f_xN_bfu(bfu(k).data) - [1-mu+r_lagrange;0;0;0];
J = f_dgdu(bfu(k).data);

c = c_0 + h*2*(rand-0.5);

cvx_begin
    variable w(m*N)
    minimize sum(abs(bfu(k).data-w))
    % minimize sum(abs(w))
    subject to
        J*w==c*g
        % abs(bfu(k).data-((1-beta)*w+beta*mmt))<=baru% Momentum
        % abs(bfu(k).data-w)<=10
        % abs(u-(beta*mmt + (1-beta)*w))<=0.15
        % norm(w,1)<= alpha_k
cvx_end